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' The unconventional Josephson coupling in a ferromagnetic weak link between d-wave super- 

conductors is studied theoretically. For strong ferromagnetic barrier influence, the unconventional 
coupling, with ground state phase difference across the link < (jigs < tt, is obtained at small crystal 
^ ' misorientation of the superconducting electrodes, in contrast to the case of normal metal barrier, 

where it appears at large misorientations. In both cases, with decreasing temperature there is an 
increasing range of misorientations, where <^gs varies continuously between and tt. When the weak 
' link is a part of a superconducting ring, this is accompanied by the flow of spontaneous supercurrent, 

^ZT^, of intensity which depends (for a given misorientation) on the reduced inductance I = 27rL/c (T) /"l?o , 

^ ' and is non-zero only for I greater than a critical value. For 1^1, another consequence of the 

unconventional coupling is the anomalous quantization of the magnetic flux. 
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O ■ I. INTRODUCTION 

Magnetic scattering effects and d-wave pairing in superconductors may have similar manifestations, such as 7r-phase 
states and spontaneous currents in Josephson devicesE 

The phase of the superconducting order parameter as a function of the momentum-space direction is different 
for d-wave and s-wave roairing states, since the latter has a single phase, whereas the d-wave state exhibits jumps 
of TT at the (110) lincs.B The Josephson current through a junction depends on the phase difference (f) between the 
, superconductors on the either side, / ~ sin <f>, and the negative sign to the prefactor of / can be regarded as an 
' additional phase shift of tt. Superconducting rings containing an odd number of such tt shifts have the spontaneous 
0^ ' magnetization ±$0/2, and quantized flux of (n -I- l/2)$07 where ri,is an integer, providing that the product of the 
On ' critical current and the self-inductance of the ring is IcL 3> $o/27r.Q 

Beside d-wave pairing, two alternative mechanisms for 7r-phase shifts at a conventional iunction have been proposed. 
One is the spin-flip scattering in the junction barrier containing paramagnetic impurities,u and the other is the indirect 
Ch ' tunneling through a localized state, in which correlation effects produce a negative Josephson coupling between two 
, superconducting grains .0 

(— I ' A problem related to the first mechanism is the possibility of 7r-coupling in aupprconductor/ferromagnetic metal 
Q \ (S/F) weak links and multilayers, due to the presence of the exchange field in F.ETlI The characteristic oscillations of 
O ' the critical tempecature Tc of multilayers as a function of the ferromagnetijcJayer thickness dp , predicted theoretically 
^ \ by Radovic et alE were observed in Nb/Gd multilayers Jay Strunk et al.EHl and by Jiang et al.Eil and in Nb/CuMn 
i (superconductor/ spinglass multilayers) by Mercaldo et al.E3 While the oscillatory Tc behavior was interpreted in terms 
/\ ' of TT-coupling in Refs. |ll| and [ij, in Ref. [l^ it was attributed to the change from paramagnetic to ferromagnetic state 
with i|*icreasing dp. Recently, double minimum Tc oscillations were observed on Nb/Co and V/Co multilayers by Obi 
et al.,Ej and the second minimum was found to be consistent with the appearance of the 7r-phase. 

In connection with the second mechanism, it was predicted that the anomalous flux quantization may appear 
with 50% probability in rings made from either disordered superconductors, or granular superconductors doped with 
paramagnetic impurities. The latter possibility was tested recently on unifornaJvIo rings doped with Fe, but the 
absence of the effect of paramagnetic impurities on flux quantization was found.llJ This may be due to the fact that 
the phase shifts induced by the impurities do not survive averaging over different electron paths around the ring. 
It was suggested that a loop with nanometer size point contact, doped with magnetic impurities, might satisfy the 
requirement that electrons circling the loop are subject to the same exchange potential. 

The Josephson effect in weak links with ferromagnetic metal barrier can be the direct probe of unconventional 
coupling. Tunneling properties-of s-wave junctions with fecromagnetic insulator barrier have been studied thepxetically 
by Kuplevakhskii and Farko,E3 De Weert and Arnold,E3 and more recently, by Tanaka and Kashiwaya.EZi In the 
latter work the possibility of tt— shift was obtained. Tanaka and Kashiwaya studied separately the Josephson, effect 
in anisotropic superconductor junctions with any symmetries and with a nonmagnetic insulating barrier.Ej This 
reference provides a general result from which several existing theories can be derived, and contains an extensive list 
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of xelfiyant references. Short weak links and d-wave junctions with normal metal barrier were studied by Barash et 
al.Eju These authors find an essentially nonharmonic current-phase relation I{(f>) at low temperature, and a flow of 
spontaneous supercurrents in superconducting ring interrupted by vr— junctions. Recently, Fogelstrom et al. have 
considered pinhole junctions in d-wave superconductors and have showia— that at low temperature the ground state 
phase difference at the junction may vary continuously between and TTr3 

In the present work, we show how new possibilities for unconventional coupling in S/F/S Josephson devices result 
from the combined effects of exchange field and d-wave pairing. We consider a Josephson weak link in the clean limit, 
with thin and short ferromagnetic (or normal) metal barrier. The interfaces between the barrier and the supercon- 
ducting electrodes are assumed fully r-transparent. Two-dimensional (2D) superconductivity and (1^2 _y2 symmetry of 
the order parameter are considered. c2l Assuming the barrier perpendicular to the a — b planes of superconducting 
electrodes, the infiuence of their relative orientation is studied. For comparison some results for 31? s-wave case 
are included. In Section II, we give a brief overview of the quasiclassical theory of superconductivity in presence of 
exchange field and for anisotropic pairing interaction, which we apply in the following. In Section III, we calculate the 
quasiclassical Green's functions for the weak link to obtain the supcrcurrent as a function of the phase difference at 
the contact. When the weak link is a part of a superconducting ring, the ground state phase difference is calculated 
and the conditions for the flow of spontaneous supercurrent and anomalous flux quantization are obtained. Section 
IV contains a discussion of the results and a brief conclusion. 

II. QUASICLASSICAL EQUATIONS 

A microscopic approach powerful enoug|h-,to deal with superconductivity in restricted geometries is the Eilenberger 
quasiclassical theory of superconductivityE3 To calculate the critical current of d-wave S/F/S contact, we generalize 
the methpd developed for the superconductor/normal metal/superconductor CS/N/S) contact by Swidzinskii and 
LikharevjHJEj applied to S/F/S contact for s-wave pairing case by Buzdin et al.El and by Demler et al.El 

The basic equations of the Eilenberger quasiclassical theory of superconductivity, in the presence of an exchange 
energy h, are 



[iw„-r3 - A + ha^i, g] + ivo • V.g = , 



(2.1) 



with normalization condition, 
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(2.2) 



where LUn = '!TT(2n+ 1) are the Matsubara frequencies {h — — 1), vq is the Fermi velocity, and g — g(vo,R;w„) 
is the quasiclassical Gorkov's Green's function integrated over energy. Explicitly, 
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the gap matrix is 



iayA(vo,R) 
i(Tj;A*(-vo,R) 

and the self-consistency equation is 

A(vo,R)^^ivr^/y(vo,v-) ^^^^^°'^'^"^y^^^^°'^''^"^ 



(2.4) 



(2.5) 



where T^(vo,Vq) is the pairing interaction. The supercurrent density is given by 

j(R)=-2ze7riVr^/vo^^^' 



(2.6) 



where (• • •) denotes the angular averaging over the Fermi surface, TV being the density of states at the Fermi surface. 



In the following, we use the notation /n = /, = /+, g| = g for the set of Green's function for one spin direction 
(down with respect to the exchange field orientation). They satisfy the scalar equations 

2{wn + ih)f + voV/ = 25A, (2.7) 

2{w,, + iK)f+ - voV/+ = 25A*, (2.8) 

voV.g- A*/-/+A. (2.9) 
For the opposite spin direction, the corresponding set of Green's fimctions is obtained by changing h —h. 

III. SOLUTIONS 

We solve the quasiclassical equations for a d-wave S/F/S junction, where S is an anisotropic superconductor with 
6,^2 _y2 symmetry, and F is a monodomain ferromagnetic metal with constant exchange energy h. 

We assume both S and F metals clean, with same dispersion relations and with same Fermi velocity vq. Electron 
scattering on impurities in S can be neglected if ^ 3> ^0, and in F if /i 3> wo/^j where I is the electron mean free path 
and ^0 the BCS superconducting coherence length. 

The thin and short barrier, of thickness d, is assumed perpendicular to the a axis in the a ~ b plane of the Ihs 
monocrystal Sl, which may be misoriented with respect to the rhs one, Sr, their d axes making an angle 9 (Fig. 
1). For anisotropic pairing, the pair potential and the shape of quasiparticle spectra depend on the misorientation. 
For d^2_y2 symmetry,E3 the pairing interaction and the pair potential in Sl are V{\'o,v'q) oc cos 293 cos 2(p' and 
A(vo) cx cos2(/3, respectively, where ip is the angle the quasiparticle momentum makes with the d axis. Similarly, in 
Sr, A(vo) oc cos2((y9 — 9). Assuming 2D nature of HTS with d-wave pairing, we take cylindrical Fermi surface for 
both metals, whereas f or 3 D s-wave case spherical Fermi surface is taken. 

To solve Eqs. (p7z|)-( |2^ ), we take in F, where h = const., A = 0, but f ^ due to proximity of S. In S, where 
h = 0, we take the order parameter on two sides of the barrier in the form Al,r = AL,R((/?)e^"^/^, where (f> is the 
phase difference at the contact, Ai^{(p) = AoCOs2(/3, Ar((p) = Aocos2((^ — 9) for d-wave (A((p) = Ao for s-wave) 
pairing Ej Assuming d/^o < Ij we neglect the spatial variation of A, which in principle could be obtained from the 
self-consistency equation (|2.5|). 

Far from the barrier, / and g approach their respective bulk values 

(/l,r) = (./i1;r)* = Al,r/c:;„l,r, (3.1) 

and 

{gL.,R) = UJ„/0JnL,R- (3.2) 

where 



'^nL.R = ^u;l + \A\l^. (3.3) 
Choosing the x-direction along cl, we look for the solution at the rhs, x > d/2, 

/r = (/r) + /Re-""", 5R = (5R)+5Re-""". (3.4) 



Similarly, at Ihs, x < d/2. 



From Eqs. (|2/7|)-(^ we find 



/l - (/l> + /Le"^^ <7L = (.9l) + .9Le"^^ (3.5) 



/r — Sr 777; /l — 3L ; 777, (3-6) 

W„-Q;RVa;/2 LUn + aLyx/2 



with 



aL,R = = — . (3.7j 

Va; Vq cos Lp 



For F, A 0, and the solutions are 

f ^ ( -2(w„ + ih)x 

= exp 



(3.8) 



and g = const, for \x\ < d/2. Note that = /(— vq, A*). Assuming a transparent boundary between two metals, we 
use the continuity conditions for g and / at the barrier interfaces, x = ±d/2. We find the normal Green's function in 
the barrier 

^ _ (Jn _ + a;nL)[(^K - a)„R)AL((/?) - {uJn - 0Jnh)^Ri'p) exp(^x)] 
WnL <I'„L[(wn - 'I;„r)Al((^) - (c^n + a)„L)AR(i^) exp(ix)] 

where 

X = </'H (3.10) 

cos (p Vo COS 

2hd 

Z= (3.11) 

is the parameter measuring the ferromagnetic barrier influence. The temperature dependence of Aq — Ao(T) is quite 
similar to the BCS one, and can be approximated by 



and 



Ao(T) = Ao(0)tanh(l.74v/Tc/T- 1) . (3.12) 



The only difference is that in the case of d-wave pairing Ao(0)/A^'^'^(0) > 10 

From the above results it is easy to obtain the corresponding ones for d-wave junction with normal metal barrier, 
putting /i — > 0, and for s-wave junction, taking isotropic pair potential A'^ — A^ = Aq, with ferromagnetic (h ^ 0) 
or normal metal (h = 0) barrier. 



A. Supercurrents 



In S/N/S case, for d-wave pairing the ground state phase difference on the contact is zero, as for the s-wave pairing, 
if two S monocrystals have the same orientation with respect to the barrier. This may not be the case for S/F/S 
contacts, and we calculate the supercurrent as a function of the barrier exchange field intensity, as well as the function 
of the orientation. 



Taking the current direction along d^, from Eq. (2.6) we get for 2D d-wave case 



v?-\Ju^ — 1 



du. 



(3.13) 



where u = l/cos(/3 and a =i,t- Green's function g^ is given by Eq. (3J), and g^{h) — gi{—h). Eq. ( 3.13| ) gives 
the supercurrent through the barrier of the area S, I = jS as a function of (p, 0, temperature T via t — T/Ao(r), 
and of parameters measuring the influence of the F barrier Z and d — d/fo(0), where fo(0) = vo/7rAo(0). Note that 
Z dft./Ao(0). The temperature dependent normalizing current is 



lo 



2Ao(r) 



(3.14) 



where the normal resistance is given by — e^voNS, N b eing t he density of states at the Fermi surface. 
In the limit T —>■ Tc, numerical calculations show that Eq. ( 3.13| ) gives, as expected, the Josephson relation 



/ = Ir sin ( 



(3.15) 



where the sign and magnitude of L. depend on T, 9 and Z . 
For 3D case and s-wave pairing,l3 Eq. (W^) reduces to 



UJn COS ■ 



■ lUJn sm ' 



9^ - 



UJn COS ^ + iujn sin ^ 



(3.16) 



where 



(3.17) 



and Eq. (|2.6|) gives 



dii, 



(3.18) 



where the same notation as in Eq. (B.13) is used. 



I supercurrents in the Josephson junctions and weak Unks are carried by the Andreev bound states in 
I In the present case, the spectrum of these states has been calculated from the analytical continuation 



In gener 
the barrier 



of the above Green's functions QcX. 



B. Magnetic flux 

For a superconducting ring containing one Joseplispn junction, the ground state phase difference (figs at the junction 
can be obtained by minimizing the reduced energycj 



VK(0) 



Jo 



/(0')d0'+-/-l(0-0,)2, 



(3.19) 



where I{4>) is calculated from Eq. ( 3.13| ), or Eq. (3.18), W — M^/(<I>o/c/27r), <I>o is the flux quantum, and the reduced 
inductance is 



$0 



with Ic = Ic{T, 9, Z). For / ^ = 0, the ground state condition reduces to I {(t)ss) — 0, I'(<i)es) _ > - 

In the limit T Tc, where the current-phase relation is harmonic, Eq. (3.15), from Eq. ( ^.19 ) it follows 



sin ipgs + I 



.) = 0. 



(3.20) 



(3.21) 



The phase (f>c — 27r$o/ $o represents the reduced external magnetic flux through the ring. In the ground state, the 
total magnetic flux 



2tt 



(3.22) 



is related to the spontaneous supercurrent Jgg = li'f'gs) by 



Note 



that 



n < $/$o < n + 1, n = 1,2,3. . ., 



^gs v^gs 



<f>n 

d) ^ $ llT 

ZTT 

now depends 
2TTn. 



on 



and 



(3.23) 
when 



IV. RESULTS AND CONCLUSION 



To illustrate the consequences of unconventional coupling we present the results of numerical calculations for typical 
cases, low temperature, t ~ 0.05 {T/Tc ~ 0.1) and high temperature, t ^ 5 {T/Tc ~ 0.9), for strong influence of the 
ferromagnetic barrier (Z — 3) and for the normal metal barrier (Z = 0). 

For d-wave pairing, /(0) dependence is much more complex than in s-wave pairing case, even in the absence of 
exchange field. Besides deformations of sinusoidal curves at low T, the shape and sign of I{4>) depend very much on 
the misorientation angle 9, both for normal and ferromagnetic metal barriers. The low T deformations of I{4>) are 
reminiscence of discontinuities at T = 0. At zero temperature, for Z = the current turns out to be discontinuous at 
(f> — TT for 6* = 0, as in the s-wave pairing case.Ej For 9 = n/2, the current jump is at = (or 2tt), and for 9 = 7r/4 



smaller jumps appear both at ^ = and tt, similarly to the results of Ref. 19, With increasing temperature, I{(f>) 
becomes less deformed tending to sinusoidal variations. Similar conclusions hold for Z = 3. The dependence I{9) 
persists in the presence of exchange field, as a strong evidence of d-wave pairing. This is shown in Fig. 2 for the 
characteristic values of Z, t and 9. 

In the absence of external magnetic field, for = the ground state phase difference (pgs at the S/F/S weak link 
is always zero or tt, both for 3D s-wave and for 2D d-wave pairing. For weak influence of the ferromagnetic barrier, 
^ ^ Ij 0gs = and for stronger influence, 1 ^ Z < 4, there is a tt— shift,|_iigs — tt (Fig. 3). Larger values of Z, for 
which 0gs — again, would correspond to the decoupling of S electrodes.EEl For 9 — 7r/2, (/)gs = — > 0gs — n and 
vice versa. These results are temperature independent. For s-wave tunnel junctions with fieiromagnctic insulating 
barrier it is also found that 0gs changes from to tt as the exchange interaction is enhanced.tZI In weak links, for 2D 
d-wave pairing we find that 4>gs varies between and tt, depending on 9 and Z. This variation is monotonous at low 
temperature, whereas at high temperature we find step function-like jumps at 9 = 7r/4. In the vicinity of 9 = 7r/4, 
depending on Z and 9 the transition from some intermediate value, < (pgs < tt, to (/)gs = or to fpis-F^ tt may occur 
with the change of temperature, similar to the to tt transition found for d-wave tunnel junction,E2l see Figs. 4(a) 
and 4(|C)|. For d-wave pinhole junctions a range of crystal orientations, where 0gs varies from to vr, is also found at 
IowTEI 

For superconducting rings with sufficiently large reduced inductance Z, interrupted by the S/F/S link, the uncon- 
ventional coupling is also found, with nonzero 0gs < tt, which rapidly decreases with increasing l^^, Figs. 4(a) and 
4(c). We note that at high temperature, t > 0.5, (pgs = for l^^ > 1. 

The appearance of spontaneous current IgsiO) in superconducting rings for external magnetic flux $e = is 
illustrated in Figs. 4(b) and 4(d) for several values of the reduced inductance. The flow of spontaneous supercurrent 
is the consequence of the misorientation effect, or of the exchange field infiuence, whenever (j)gs ^ 0. This is a 
generalization of the spin-flip induced effect predictedj-by Bulaevskii et al. for conventional superconductor rings with 
junction barrier doped with paramagnetic impurities.u The highest values of spontaneous currents, corresponding to 
9 = tt/2 for Z = 0, and to 9 = for Z = 3, strongly depend on the reduced inductance. Maximum values, equal to 
the critical currents, correspond to l^^ ~ 1 at low T, and to < 1 at high T (Fig. 5). It is important to point 
out that the spontaneous currents appear only for sufficiently large reduced inductance, greater than some critical 
value Ic, in contrast to the result of Barash et al. for S/N/S junctions at T = O.Ej At high T, t > 0.5, where tha 
current-phase relation is harmonic, we obtain an universal curve below l^^ — 1, as in the case of tunnel junctions.u 
At low T, the spontaneous currents flow is obtained in a wider range of the shape of the curves and Ic depending 
on the barrier influence Z and on the type of pairing symmetry. 

In the presence of external magnetic flux, $e 0, the flow of spontaneous current leads to the anomalous flux 
quantization in the ring, both for S/N/S and S/F/S junctions. In the first case, Z — 0, this occurs in the vicinity 
oi 9 = tt/2 and in the second case, Z = 3, in the vicinity of 6* = 0. The effect of half magnetic flux quantization, 
$/$o = 1/2,3/2, . . ., pronounced for small, becomes smeared out and eventually lost for larger (Fig. 6). 

Since oc l//c(T) rapidly increases with increasing temperature, the unconventional coupling effects, such as the 
spontaneous supercurrent flow and the anomalous flux quantization exist only below some temperature characteristic 
for the given device. 

In conclusion, the spontaneous supercurrent flow in a HTS ring with ferromagnetic weak link provides new pos- 
sibilities for experimental observation of the unconventional Josephson coupling, < 0gs < tt. Due to the presence 
of the exchange field in the barrier, it may appear, manifested by the anomalous flux quantization, without the 
misorientation of the crystals in two superconducting electrodes. However, the exchange field effect do not mask the 
difference between s and d-wave pairing, since the coupling depends on the misorientation. We emphasize that the 
unconventional coupling effects in the weak links are much more diverse than in the tunnel junctions, which display 
in the ground state only or tt phase difference. 
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FIG. 1. Schematic illustration of the weak link. The fourfold symmetry of the gap functions is indicated. 



FIG. 2. Reduced current I/Io as a function of phase difference (j>, for S/N/S, Z = Q and S/F/S, Z = weak links, at low 
temperature t = 0.05, and at high temperature t = 5. The reduced thickness is d = 0.25. 2D d-wave pairing for three values of 
the misorientation angle (solid curves) , and 3D s-wave pairing (dashed curves) . 



FIG. 3. Ground state phase difference across the link 0gs as a function of the ferromagnetic barrier influence Z, for = and 
d = 0.25. 2D d-wave pairing (solid lines), and 3D s-wave pairing (dashed lines). These results are temperature independent. 

FIG. 4. Ground state phase difference <f>gs and reduced spontaneous current /gs//c as functions of misorientation angle 6 
for a d-wave superconducting ring interrupted by S/N/S, Z = and S/F/S, Z = 3 weak liuks, with d = 0.25 and for three 
values of the reduced inverse inductance, = (open ring), 0.1 and 1. Low temperature, t = 0.05 (solid curves) and high 
temperature, t = 5 (dashed curves). 



FIG. 5. Highest spontaneous supercurrent |7gB|//c as a function of l~ in a superconducting ring interrupted by a weak link. 
Reduced barrier thickness d = 0.25. Low temperature, t = 0.05 (solid curves): {!) Z = 0, 9 = it/2 and (2) Z = 3, 6 = for 
2D d— wave pairing, and (3) Z = 3 for 3D s-wave pairing. High temperature, t = 5 (dashed curve): universal behavior for all 
above cases . 

FIG. 6. Magnetic flux $ through a d-wave superconducting ring as a function of the misorientation angle 9 ior t = 0.05, 
d = 0.25, for two values of and four values of the external magnetic flux $e. S/N/S, Z = 0, and S/F/S, Z = 3 weak links. 
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